Quintessence with quadratic coupling to dark matter 



Christian G. Bonnier, 1 ^ Gabriela Caldera-Cabral, 2 '0 Nyein Chan, 1 '© Ruth Lazkoz, 3 'EI and Roy Maartens 2 '© 

1 Department of Mathematics and Institute of Origins, 
University College London, London WC1E 6BT, UK 
2 Institute of Cosmology & Gravitation, University of Portsmouth, Portsmouth POl 3FX, UK 
J Fisika Teorikoa, Euskal Herriko Unibertsitatea, 48080 Bilbao, Spain 
(Dated: November 16, 2009) 

We introduce a new form of coupling between dark energy and dark matter that is quadratic 
in their energy densities. Then we investigate the background dynamics when dark energy is in 
the form of exponential quintessence. The three types of quadratic coupling all admit late-time 
accelerating critical points, but these are not scaling solutions. We also show that two types of 
coupling allow for a suitable matter era at early times and acceleration at late times, while the third 
type of coupling does not admit a suitable matter era. 
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Cosmological observations strongly suggest that the expansion rate of the universe is accelerating and that matter 
in the universe is dominated by non-baryonic cold dark matter (see e.g. [l|). However, what exactly causes this 
acceleration is not well understood, and one of the main challenges of modern cosmology is to understand the nature 
of this mysterious dark energy. The existence of some form of dark matter is long known, as implied by the flattened 
galactic rotation curves observed by Zwicky as early as 1933. Several experiments have been carried out (see, for 
example Q) in search of candidate dark matter particles. The fact that dark matter only interacts weakly with 
standard matter means that it is difficult to detect such particles directly. Neither dark energy nor dark matter have 
been detected directly. Only the total dark sector energy-momentum tensor is known from its combined gravitational 
effect. In order to separate the two components, we have to assume a model for them. It is possible that these 
components interact with each other, while not being coupled to standard model particles. Such a possibility can 
lead to new approaches to the coincidence problem ( "how do dark matter and dark energy attain the same order-of- 
magnitude value at the right time to allow for the observed large-scale structure?" ) . It can also produce interesting 
new features in large-scale structure, such as a large-scale gravitational bias [|[ and a violation of the weak equivalence 
■ principle by dark matter on cosmological scales [J] . 

In this paper we study a class of cosmological models with interactions in the dark sector. Various models of the 
coupling between dark energy and dark matter have been proposed and investigated (see, e.g. [S [S 0, H, H, EE EE 
EE [H, EI EE EE EE El, ESM EH ) . We consider only the background dj r namics - for cosmological perturbations of 
coupled dark energy models, see e.g. 0, M, US US Hi US US HE US HI, M, US EE US EH HE IM HE EE US EE EE 

E3,E|. 

There is no fundamental theory that selects a specific coupling in the dark sector, and therefore any coupling 
model will necessarily be phenomenological, although some models will have more physical justification than others. 
Here we analyse the background dynamics for a new model of coupling. This model improves the one previously 
introduced in [TEl . l2lj . which was motivated by simple models of inflaton decay during reheating and of curvaton 
decay to radiation. 

The background description of a coupled model with quintessence dark energy density p v and dark matter density 
p c is given by the energy balance equations 

p c = -3Hp c + Qc, (1) 
pip = —3H(1 + w v )p^ + Q v , Q ip = -Q C :=Q. (2) 
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Here Qa is the rate of energy transfer to species A. It follows that 



, . > r • f dark matter — > dark energy 

O energy transfer is < , , , , , , (3) 

< 1 dark energy — > dark matter v ' 



The dark energy equation of state parameter is 



_ Py _ l^ 2 ~ V{<P) 



(4) 
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The modified Klein-Gordon equation follows from Eq. ((2]) as 



dV Q 

For quintessence with an exponential potential, 

V(tp) = V exp (-nXtp) , k 2 := 8nG , (6) 

where A is a dimensionless parameter and Vq > 0. We neglect the radiation and therefore the evolution equations are 

p b = -3Hp b , (7) 
k 2 

H = -—[p c + Pb + tp 2 ] , (8) 
where baryons pb are not coupled to the dark sector. The Friedman constraint is 

h2 = Y (i ° c + pb + Pv) ■ (9) 

We define effective equation of state parameters for the dark components which describe the equivalent uncoupled 
model in the background: p c + 3H(1 + w CtC ff)p c = 0, p v + 3H(1 + W Vte s)p<p = 0. By Eqs. (JTJ) and ©, 

W c ,cS = TTrr ) w V,off = W V - — . (10) 

0/1 p c ■in p v 



II. NEW MODEL OF DARK SECTOR COUPLING 

In [2l| a model of the form 

Q = r Pc , (11) 

was introduced, with T constant. The motivation for this form of interaction is that, for T > 0, the same Q is used 
for simple models of: (1) the decay of an inflaton field to radiation during reheating [44j, (2) the decay of dark matter 
into radiation [45||, (3) the decay of a curvaton field into radiation |46[ [. (4) the decay of super-heavy dark matter 
particles into a scalar field [47j |. 

We consider this coupling to be better motivated than alternatives of the form Q = aHp c - which are designed for 
mathematical simplicity, since they lead to the same number of dimensions of the phase space (two) as the uncoupled 
case. The coupling in Eq. (fTTj) is not designed for mathematical simplicity, but is chosen as a physically simple 
form of decay law. It leads to a three-dimensional phase space. This new phase can be compactified [15| . as in 
the two-dimensional case, but great care is required in analysing the stability properties of the resulting dynamical 
system. The stability matrix contains singular eigenvalues as one approaches the critical points. In [l5[ we developed 
the required machinery to overcome these problems and were able to present a complete phase space analysis. Our 
techniques are readily applicable to more general couplings. 

Simple decay laws of the form in Eq. (fTTjl fail to reflect the feature that interactions are typically determined by 
both energy densities. We therefore consider the natural first extension Eq. (fTl]) to a quadratic form 



Q = Ap 2 v + Bp 2 c + Cp c p v , 



(12) 
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where A, B and C are coupling constants. We define dimensionless coupling constants as 

a = AH , (3 = BHq , j = CH Q . (13) 
The Friedman constraint (J9j) in dimensionless form becomes 

n c + n„ = i, fi: =|^> ( 14 ) 

where we neglect the baryons, and the total equation of state parameter is given by 

p %ot + 3H(1 + to t ot)ptot = , (15) 
totot : = = — — — = w v Qm ■ (16) 

Ptot P V + Pc 

The condition for acceleration is w to t < —1/3. A phantom field with w v < ~ 1 violates the dominant energy condition, 
p > \p\. We therefore assume that w v > — 1, thereby excluding phantom models with negative kinetic energy. 
We introduce the dimensionless variables x, y, as in the uncoupled case (48|, where 

Then y > because of the positivity of the potential energy, and Eq. (fl4|) implies that 

< n v = x 2 +y 2 < 1 . (18) 
In the new variables, the equation of state parameters are 

x 2 2 

w v = — 9 i ^tot = x 2 - y 2 . (19) 
x 1 + y 2 - 

The Hubble evolution equation may be written as 



IP ~ ~2 



2 =-^(l + x 2 -y 2 ). (20) 



As already indicated, it turns out that the resulting evolution equations do not allow a two-dimensional representation 
of this model, since we cannot eliminate H from the energy balance equations ([T]) and using only the variables 
x(N),y(N), where we use N — log(a) as the independent variable. Equation (|20|) must therefore be incorporated into 
the dynamical system. We do this via a new variable z, chosen so as to maintain compactness of the phase space 

UU (2D 



H + H< 



Thus < z < 1, and the compactified phase space now corresponds to a half-cylinder of unit height and radius [l5j. 
The top of this half-cylinder is defined by z = 1 and it turns out that the equations become singular as z — > 1. 
Therefore, care is required in order to analyse the resulting dynamical system. 

It is evident that the quadratic and higher-order couplings can be treated in a similar fashion, for instance one 
could consider a coupling of the form Q = pp^p c . The most general model is of the form 

Q = $>m»prps, ( 22 ) 

where m, n are non-negative integers. The matrix q mn is arbitrary except for the condition q^Q = 0. Note that q mn 
has no a priori symmetry properties and is not necessarily a square matrix. The linear model in Eq. pip was extended 
to the most general linear model, Q = —(T c p c +T t j,p < j > ) in [49J. The linear and quadratic models lead to 

(ff)mn= (° T Q ) >IJt ^ ! ° C " I ' (2 - !) 
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III. DYNAMICAL ANALYSIS 



In this section we analyse the three particular cases when two of the interaction terms are equal to zero. Two of 
these models allow for a standard matter era, but the model with A = C = does not allow it. Then we combine the 
models A and C and analyse the composite model. 



A. Model A: coupling Q = -§^P V 
The system of autonomous differential equations is 

/ n xV^ 2 3 . 2 2 N 3(1 ~ z)(x 2 + y 2 ) 2 
x' = -3x + A-y y 2 + -x(l + x 2 - y 2 ) + a — 



2xz 



I \ 3 . 2 2 \ 

V = ~ X ~ X V + +x ~ V ' ' 



z' = -z(l - z)(l + x 2 



v 2 ). 



(24) 

(25) 
(26) 



The critical points, defined by x' — 0, y' = and z' = 0, and the eigenvalues of the stability matrix are given in Table[U 
In Table ITTl we characterize the critical points and give the effective equation of state for the late-time attractor. 



Point 


£* 


y* 


z* 


Eigenvalues 


A 











3 3 3 
2 ' 2 ' 2 


D 








1 


3 3 3 
2 ' 2 ' 2 


E± 


±1 





1 


3 , 3 , 3 -p y^^" ^ 


F 


ps l 

V 2 A 


l~3 1 
V 2 A 


1 


-f ,-^(A±V24-7A2) 


G 


A 

76 




1 


4r 7 — 3 + , — 3 + A 



TABLE I: Critical points and associated eigenvalues for coupling model A. 



Point 


Stable? 




WT 


Acceleration? 


Existence 


A 


Saddle node 








No 


VA, a 


D 


Saddle node 








No 


VA, a 


E± 


Saddle node 


1 


1 


No 


VA, a 


F 


Stable focus for A 2 > ^ 
Stable node for 3 < A 2 < ^ 


3 
1? 





No 


A 2 > 3 


G 


Saddle node for A 2 > 3 
Stable node for A 2 < 3 


1 


- — 1 

3 1 


A 2 < 2 


A 2 < 6 



TABLE II: The properties of the critical points for model A. 



This model depicts an evolution of the universe in good agreement with the observations for certain values of 
the parameters A, a. Saddle point A corresponds to the standard matter dominated universe with a(t) oc i 2 / 3 ; its 
instability allows for the existence of trajectories escaping from it and ending at an attractor, which exists for adequate 
values of the parameters. However, the attractor (or late time stage of the universe) will only represent an accelerated 
scenario for a flat enough potential, specifically when A 2 < 2. In this case the attractor is completely dark energy 
dominated, point G. If the potential is not flat enough, the attractor, point F, is a scaling solution, in which the 
fraction of dark energy is the dominant one. 
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One way to understand this behaviour is by looking at the relative strength of the coupling, /, in the CDM balance 
equation, 

/-■g-. (27) 

In the matter era, H 2 — p c /3, we have 

f-pl/pT, (28) 



which is decreasing into the past. Thus the coupling is weaker in the past and this allows a near-standard matter 
era. It is important to note that the dynamical behaviour of the standard matter point A and the accelerated critical 
solution G does not depend on the sign of the coupling parameter a. 

B. Model B: coupling Q = -j^pl 

The autonomous system is 

/ n 2 3 . 2 2 n „3(1- z)(l-x 2 -y 2 ) 2 
x 1 = -3s + y 2 + -x(l + x 2 - y 2 ) + ^ , (29) 

y' = -X^x y + ^y(l + x 2 -y 2 ) 1 (30) 

z' = \z{l-z)(l + x 2 -y 2 ). (31) 

The critical points and their stability properties are summarized in Tables IIIII and IIVI respectively. 

This model does not have a unstable standard matter solution, so it cannot depict an evolution from an early 
dustdike scenario. However, there is a solution in which the dark energy can mimic such behaviour for A 2 = 3. This 
is saddle point C, but the A value allowing for that point to act like matter prevents the existence of an accelerated 
solution at late times. Regarding the late-time attractor we meet again the same situation as before, a scaling but 
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Point 








Eigenvalues 


Wtot 


R _l 




n 

u 


n 

u 


^ 3 3 =n . /- X 
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Kj 


A 

Ve> 


v 1 -- 


U 


\" A 2 q \2 O 

2 ' 2 J ' A ° 


— - 1 


t\ 
U 


u 


u 


1 


3 3 3 
2 ' 2 ' 2 


u 


F j_ 


±\ 


n 
u 


\ 


_Q Q Q T- , / 2 \ 
G , O , O -(- A / 2 


1 


F 


fz 1 
y 2 a 


[Jl 

V 2 A 


1 


-f ,-^(A±V24-7A2) 





G 


A 

%/6 




1 


4j- , 3 + 4|- , 3 + A 


V- 1 



TABLE III: Critical points and associated eigenvalues for coupling model B. 



Point 


Stable? 


"</> 


Wt 


Acceleration? 


Existence 


B+ 


Saddle node for A > v6 
Unstable node for A < v6 


1 


1 


No 


VA, /3 


B_ 


Unstable node for A > — \/6 
Saddle node for A < — v6 


1 


1 


No 


VA, p 


C 


Saddle node 


1 


h 1 

3 1 


A 2 < 2 


A 2 < 6 


D 


Saddle node 








No 


VA, p 


E± 


Saddle node 


1 


1 


No 


VA, p 


F 


Stable focus for A 2 > 21 
Stable node for 3 < A 2 < 24 


3 





No 


A 2 > 3 


G 


Saddle node for A 2 > 3 
Stable node for A 2 < 3 


1 


- — 1 

3 1 


A 2 < 2 


A 2 < 6 



TABLE IV: The properties of the critical points for model B. 



dark energy dominated non accelerated solutions for excessively shallow potentials, and an accelerated completely 
dark energy dominated scenario in the opposite case. The relative strength of the coupling Eq. (f2"T)) for this model is 

f~H, (32) 

and is increasing to the past. Since the coupling gets stronger at early times we cannot get a standard matter era. 
The direction of the energy exchange does not affect this conclusion. 



In this case, 



C. Model C: coupling Q = j^p c p v 



, , ,V6 2 3 2 2 i{l-z)(l-x 2 -y 2 ){x 2 + y 2 ) 

x = -?,x + \— y +-x{l + x -y + 7 , (33) 

2 2 2xz 

l/ = -X^-xy+^y(l + x 2 -y 2 ), (34) 
z' = ^z(l~z)(l + x 2 -y 2 ), (35) 



and we summarize the critical points and stability in Tables [V] and I VII 

This model has interesting properties. First of all, a standard matter represented by the unstable point A exists. In 
particular when 7 > 0, the instability of the point is generic, as it is occurs in all directions; whereas for 7 < the point 



FIG. 2: Phase-space trajectories for model £>. The left plot shows the stable node G, with A = 1.2 and f3 = 10 3 . The right 
plot shows the stable focus F with A = 2.3 and (3 — 10~ 3 . 



Point 


a;* 


y* 


2* 


Eigenvalues 


Wtot 


A 











| , | ,sgn(7)oo 





B± 


±1 








3 ' 3= F y§A,-sgn(7)oo 


1 


C 


A 







^- - 3,-sgn(7)oo 


1 


D 








1 


3 3 3 
2 ' 2 ' 2 





E± 


±1 





1 


3,3,3 — (— \J~^ A 


1 


F 


fz 1 

V 2 A 


/¥i 

V 2 A 


1 


-f ,-^(A±V24-7A^) 





G 


A 




1 


-^,-3 + ^,-3 + A 2 


\2 

- — 1 

2 x 



TABLE V: Critical points and associated eigenvalues for coupling model C. 

becomes a saddle. Secondly, an accelerated attractor exists for adequate values of the parameters. The two possible 
late-time attractors as in the other two cases are found again, and the same conditions as in those cases operate in 
connection with the kinematical features of the scenario they represent, if A is not small enough, acceleration will not 
be possible and in addition dark energy will not dominate completely. 

The coupling strength for this model, Eq. ([??]) is decreasing into the past as 

/ - P* 2 /P 1 J 2 ■ (36) 

D. Superposition of couplings 

When we combine the different coupling models together, we expect only those critical points to be present which 
are critical points of each individual model. Since we want to describe the evolution of a universe that includes a 
standard matter era and evolves towards a stable accelerating solution, we must choose the combination of models 
A and C, since those are the ones that allow for a standard matter era. Therefore we consider the composite model 
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Point 


Stable? 


Q4, 


WT 


Acceleration? 


Existence 


A 


Unstable node for 7 > 
Saddle node for 7 < 








No 


VA, 7 


B+ 


Unstable node forA < and 7 < 
Saddle otherwise 


1 


1 


No 


VA, 7 


B„ 


Unstable node for A > —y/E and 7 < 
Saddle otherwise 


1 


1 


No 


VA, 7 


C 


Saddle node 


1 


£-1 


A 2 < 2 


A 2 < 6 


D 


Saddle node 








No 


VA, fj 


E± 


Saddle node 


1 


1 


No 


VA, fj 


F 


Stable focus for A 2 > f 
Stable node for 3 < A 2 < ^ 


I? 





No 


A 2 > 3 


G 


Saddle node for A 2 > 3 
Stable node for A 2 < 3 


1 




A 2 < 2 


A 2 < 6 



TABLE VI: The properties of the critical points for model C. 

E- 




FIG. 3: Phase-space trajectories for model C. The left plot shows the stable node G, with A = 1.2 and 7 = 10 . The right 
plot shows the stable focus F with A = 2.3 and 7 = 10 -3 . 



defined by the coupling 



Q = IT Pip + -FTPcPv ■ 
-Ho -Ho 



(37) 



We note that the two couplings are decoupled in the sense that there are no cross coupling terms in the dynamical 
system. 



9 



The resulting system of autonomous differential equations reads 

o 2 3 2 2 3(l-z)(x 2 +y 2 ) 2 

x = -3x + A— y + -x(l + x -y ) + a 

2 2 2xz 

, v6 3 2 2\ 

V = -^—xy + -y{l + x -y), 

z' = ^z(l-z)(l + x 2 -y 2 ). 
The critical points and their stability are listed in Tables I VIII and IVIIII 



3(1 - z)(l - x 2 - y 2 ){x 2 + y 2 ) 
2xz 



(38) 

(39) 
(40) 



Point 


x* 


y* 


Za. 


Eigenvalues 




A 











sgn( 7 )oo , | , | 





D 








1 


3 3 3 
2 > 2 ' 2 





E± 


±1 





1 


-3,3,3 T ^/|A 


1 


F 


i fz 
Ay 2 


i fz 

A V 2 


1 


-f ,-^(A±V24-7A2) 





G 


A 




1 




^ - 1 

2 x 



TABLE VII: Critical points and associated eigenvalues for the superposition of couplings for model A and model C. 



Point 


Stable? 




Wt 


Acceleration? 


Existence 


A 


Unstable node for 7 > 
Saddle node for 7 < 








No 


VA, a, /3, 7 


D 


Saddle node 








No 


VA, a, /3, 7 


E± 


Saddle node 


1 


1 


No 


VA, a, A 7 


F 


Stable focus for \ 2 > ^■ 
Stable node for 3 < A 2 < f 


3 

A^ 





No 


A 2 > 3 


G 


Saddle node for A 2 > 3 
Stable node for A 2 < 3 


1 




A 2 < 2 


A 2 < 6 



TABLE VIII: The properties of the critical points for the superposition of couplings. 
Points A and D correspond to the standard matter era and point G is the accelerated attractor for A 2 < 2. 

IV. CONCLUSIONS 



We have made a comprehensive analysis of the background dynamics for a new class of models with quadratic dark 
sector coupling, which are a simple physically-motivated generalization of the linear coupling couplings with constant 
rates of energy transfer, given in 

[TEHllii. The tw ° 

species which interact are dark matter and quintessence with 

an exponential self-interaction potential. 

Of the three different terms in the general quadratic coupling, we found that the term Q = Bp 2 leads to a universe 
without a standard matter era, whereas the other two terms, Q — Ap^ and Q — Cpcp 2 ,, do admit a standard matter 
era and an evolution that connects this to a late-time attractor. This attractor is accelerated provide the potential is 
flat enough. The models we have analysed provide us with the following partial answers. These features are valid for 
both signs of A or C, i.e the evolution is not affected by the direction of the energy transfer. But in the C > case 
the instability of the matter era is more generic; so there is in a way more room for a transition from the matter era 
to the accelerated attractor. In other words, in theses case there are less restrictions on the initial conditions for this 
desired transition between asymptotic states to occur. 
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The critical point for late-time acceleration, G, is not a scaling solution, since 

o c « = o, n y * = i. (41) 

This is similar to the asymptotic behaviour of the standard ACDM model. Thus the quadratic models do not produce 
a constant non-zero and finite ratio r2 c */f20*, and therefore do not address the coincidence problem in this sense. The 
linear coupling Q = Tp c leads the same critical point G; this accelerated solution is an attractor when T > 0, i.e when 
dark matter is decaying into dark energy. 

The quadratic models which admit a viable background evolution can be compared to observations in order to 
constrain the parameters a and 7. This will require an investigation of the cosmological perturbations in these 
models. 
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